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Abstract: In this paper a new saturation model is presented. This model is based on the the- 
oretical solution for the generating functional, and it is quite different and not more complicated 
than the Glauber-like approach used before. The model describes the structure function F 2 of the 
proton, as well as the diffractive structure function ■ We show the difference between our model 
and the eikonal approach by calculating the multiplicity distribution, using the AGK cutting rules 
strategy. 
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1. Introduction 

The very intriguing observation of the HERA experiments is the rapid rise of the total 7*p cross 
section, with energy in the deep inelastic scattering (DIS) region. The linear evolution equations 
approach is not valid in the region of high energies, since it predicts an increase of the parton 
density, steeper than is allowed by the Froissart bound. The solution to this problem is hidden in 
the effect of the parton saturation. At high energy the density of partons increases, filling the whole 
transverse area of the target. This stage is called saturation. At small values of Bjorken-x, the 
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interaction between partons should be taken into account leading to the recombination of partons. 
This process tames the growth of the parton density. 

The main goal of this paper is to build a model which follows on from the self consistent theoret- 
ical approach, and is able to describe the matching between soft and hard processes. This approach 
is based on the generating functional method. The model incorporates the following features : 

i) Evolution of the parton density with energy 

ii) The phenomena of parton saturation. The model takes into account the effect of parton 
recombination at high density and preserves unitarity. 

Hi) This is a simple model, different from the model based on the eikonal approach. 

iv) Good description of all recent experimental data on inclusive, hard, soft, and diffractive 
processes. 

In the next section, we start with a brief discussion of the main properties of the generating 
functional [1]. We find the solution to the simplified evolution equation, assuming that the dipoles 
do not change their transverse sizes, during the interaction. 

Based on this solution, we build our model. All details and properties of the model are discussed 
in section ^| This model is quite different from the other saturation models, which are based on 
the eikonal approach [2-7] . In section |3|, we show how our model describes the HERA data using 
the recent data on the proton structure function F 2 , in a wide range of kinematics from different 
collaborations [8-16], we fit our model and find the free parameters. We compare the developed 
model, to other saturation models in section f|, where we discuss the multiplicity distribution of the 
two models, and we present the differences between them. In the next section using the parameters 
from the fit, we describe the diffractive dissociation experimental data, the charm quark structure 
function F| c , and the slopes (IF2/ d(\riQ 2 ) and dlni^/rf(ln 1/x). Finally, we summarize our results. 

2. Description of the Model 

2.1 General properties of the generating functional 

This approach is based on the concept of color dipoles [1], which in the large N c limit are correct 
degrees of freedom at high energies [17], since they diagonalize the scattering matrix. At high 
energy, every parton tends to emit more partons, which leads to the evolution of the initial parton 
density. The dipole decay produces a parton cascade, originating from the one parent dipole at 
initial high rapidity Y. We define by P n to be the probability to obtain n dipoles from the parent 
dipole, at some rapidity y. 

Using the generating functional approach, one can construct an evolution equation for the dipole 
probability P n at rapidity y [18]. The interesting property of this approach, is the manifestation 
of the saturation phenomena, and the conservation of unitarity. For the sake of simplicity, we 
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consider the generating functional approach in the toy-model [18,19], where we neglect changes to 
the dipole sizes during interaction. In this approach, we obtain the following evolution equation for 
the probability P n 

5— = -UanPn + uj (n - l)P n -i (2.1) 

oy 

where u Q is the probability for one dipole to decay into two. The generating functional, with 
our assumption degenerates into the function 

Z(y,u) = p n(y)u n , (2-2) 

n 

Equation Eq. ( |2.1|) , for the probability P n can be rewritten in terms of the generating functional 
Eq. (13), as 

9Z{y, u) dZ(y, u) 

= -u u{l - u) — (2.3) 

oy ou 

This is a Liouville equation which has the solution 

u e -<*>o (y-yo) 

Z(y -y,U) = ; ; ; r . (2.4) 

To find the interaction amplitude, we follow the procedure suggested in Ref. [20], namely, 

00 1 f) n 7 00 1 f) n 7 

N(y, b) = ImA^y, b) = - £ - t^I^i ln(b) = - 7^U=i iW) (2.5) 

^— ' n! (ou) n n\ (ou) n 

n=l v ' ra=l v ' 

Here, — 7 n (6) is the imaginary part of the amplitude for the scattering of n dipoles off the target, 
at fixed impact parameter b, at low energies. The most important assumption, is the independent 
interaction of n dipoles with the target, which is expressed by the factor 7„ = 7™ factor in Eq. (|2.5| ). 



2.2 The Model 

Using the generating functional approach, we propose a model for the interaction amplitude of the 
dipole with the target. The amplitude N(y, b) can be found from the following relation 

N(y; b,r) = 1 - Z(y, u(r)) . (2.6) 

Inserting Eq. (|2.4j) into Eq. (|2.6|), one obtains the final expression for the interaction amplitude 



N(v b r) = — { - -- = 71 ' (2 7) 

K ' ' (1 - u)e?<>(v-y°) + u 7(r)e a; o(»-fo) + (1 — 7(r)) v > 
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Eq. ( |2.7| ) describes the sum of so- called" fan" diagrams. Introducing the new variable [21] 



1R = (2.1 
1-7 



the amplitude now has the simple form 



N(y, b, r) = ^^^)f ( 7° } (2 . 9) 

and satisfies the following equation 

dN(y;b,r) dN(y; b, r) 

This means that in fact, the generating functional describes the system of non-interacting hard 
pomerons, but its interaction with the target should be renormilized using Eq.( |2.2| ). Our main 
assumption, is that we can include the dependence on the size of the dipoles incorporates it in the 
low energy amplitude 7 condition. 

1r{Vo\ b) -> 7tf(y ; b, r) (2.11) 

where r is the transverse size of the dipole. At small values of r, the amplitude of Eq. ( |2.9| ) should 
match the expression originating from the perturbative calculation of the BFKL pomeron exchange, 
namely 

Q = 7 *(yo; b, r)e"°^ = ^ r 2 a s {^)G{y ^ 2 )G{y -y ,fi 2 )S(b) (2.12) 



where G(y — y , /i 2 ) is the gluon density at some scale /i 2 , y is the rapidity defined in Eq. ( |2.12| ) as 
y = ln(l/x), and S(b) stands for the proton profile function in the form 1 

*>-^(£)*.(£l 

where K\ is the McDonald function. The evolution of the parton density is given by the DGLAP 
evolution equation, with the initial parton density at some scale Q\ 

G(x,Ql) = ± (2.14) 
where A and uj are to be determined from the data fit. The expression for the gluon density 



1 This S(b) is the Fourier transform of the dipole formula for the electromagnetic form factor of the proton. We 
can use this as the first approximation, having in mind that the real 6-dependance should be determined by the 
impact parameter dependence of the low energy amplitude for dipole - proton scattering, so-called "two gluon form 
factor" . 
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G(x, Q 2 ), is obtained from the solution to the DGLAP equation, and is given by the inverse Mellin 
transform [22] 

I rc+ico ^ 

G{y,t) = — duexp{--t + uy)- (2.15) 



where the variable t is defined as 

4iV Cl ln(Q 2 /A 2 ) , f 2n fs , . . . . . 

where Qq is the initial condition, and rif is the number of flavors. The analytical solution to 
Eq. (|2.15| ) reads as 



G(y,t) = Ae~ t+UJoy 



(2.17) 



Here, we introduce the hard scale /i 2 , which corresponds to the transverse size of the dipoles in the 
following way 2 



g 2 ^/i 2 = ^ + /i 2 (2.i8) 



where we rewrite the gluon density in terms of fi 2 , rather than Q 2 . The parameters C and //q, 
are obtained from the fit to data. Finally, the expression for the interaction amplitude takes the 
following form 

with 

n(r,b,x) = ^r 2 a s (/i 2 (r))G(ln(l/x ),/i 2 (r)) G (\n(x/x ), fi 2 (r)) 5(6) (2.20) 

This is a new form for the interaction amplitude of the dipoles with the target, which originates 
from the generating functional approach. 

2.3 Saturation Scale 

At some energy value, partons start to populate densely, and this leads to the effect of saturation. 
We define a new scale, which separates the two regions, namely, the region of low parton density, 
where we can apply perturbative methods, and the region of high parton density, where we should 
take into account the recombination effects, and non-liner corrections to the parton density. In order 
to determine the saturation scale, we demand that the packing factor of the partons Eq. (|2.21|) , at 



2 Such a form of parametrization was used in [2-7], and we introduce the same parametrization, to make the 
comparison easier between our model and other models. 



- 5 - 



some energy, be equal to 1. The packing factor is denned as 



3_ 

TlR 2 



(2.21) 



where <jq is the typical interaction cross section of the partons, which is proportional to Q 2 . N g 
corresponds to the number of partons, and 7r R 2 denotes the area of the transverse slice of the 
hadron. Putting everything together, one obtains 



1 = j^r 2 a s (fi 2 s )G(y ,^ s )G(y -y ,fi 2 s )S(b) 



Q 2 N C 



7T 

Ql( x ) = ^jf^s(fil)G(y , t^)G(y - y , fi 2 s )S(b) 



(2.22) 



where fi 2 = + /j,q. The estimation of the saturation momentum is obtained from the numerical 



solution of the Eq. ( 2.22j ) by iteration. The result obtained for the saturation scale is plotted in 



Fig. 2.1 
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Figure 2.1: Saturation scale for various values of impact parameter h(fm) as a function of Bjorken-x. 



Different graphs correspond to different values of the impact parameter b. The value of the 
saturation scale, is calculated in units of GeV 2 . 



3. Description of DIS 

We start by overviewing the main features of DIS. The deep inelastic scattering process is shown in 
Fig. |3.1| (a) and the standard variables are defined as 
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P(P) 



(a) (b) 
Figure 3.1: Kinematics of deep inelastic scattering and photon-gluon fusion diagram. 



incoming proton momentum : P 

photon's virtuality : Q 2 = -q 2 

P-q 



fraction of electron energy transferred to the proton : y 

em energy : (( 

Bj or ken- x : x 



P-k 

photon - proton system energy : (q + Pf = W 2 

Q 2 Q 2 



2p-q Q 2 + W 2 



(3.i; 



For small values of Bjorken-x, one can separate between the photon wave function, which corre- 
sponds to the creation of the quark - antiquark pair, and the interaction of the dipole with the 
target [23,24]. In the dipole picture of DIS, the total photon-proton scattering cross section from 
transverse (T) and longitudinal (L) polarized virtual photons, is given by the convolution of the pho- 
ton wave function, ^t,l, and the interaction amplitude N of the dipoles with the target Eq. ( p,19| ). 
The proton structure function F2, can be written as the sum of two contributions [24] 

F 2 (x,Q 2 ) = 2^—^ f d 2 r [ d 2 b [ dz{\y T (r,z;m f ,e f )\ 2 + \V L {r, z; m f , e f )\ 2 } N{r, b, x) 

2 47T a em J J Jo 

(3.2) 

where / stands for the flavor of the quark-antiquark pair, and rrtf and e/, denote the mass and 
electric charge of the quark with flavor /, respectively. The squared wave function of the photon is 
given by 

N e 2 a 

\y T (r,z;m f ,e f )\ 2 = ° ^ {[z 2 + (1 - z) 2 ]a 2 K 2 (ar) + m 2 f K 2 (ar)} (3.3) 



and 



\* L (r,z;m f ,e f )\ 2 = { 4 QV(1 - z) 2 K 2 (ar)} (3.4) 
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where 



a 2 = z(l - z)Q 2 + m) (3.5) 

and K%, Kq where K\, Kq are McDonald functions. The significant contribution to the inclusive 
processes, appears from the heavy quark production. In the present paper we investigate the 
influence of the charm quark contribution to the proton structure function F 2 . 



3.1 Heavy quark 

The recent data from HERA, shows that the heavy quark contribution is up to 30%-40% of the 
total value of the structure function [11,12], and can not be neglected. Thus, we take into account 
the contribution of the charm quark and confront it with recent data. The analysis does not include 
any new parameters. The proton structure function may be written, as the contribution from all 
possible participating quarks. Since the contribution of each quark flavor is proportional to its 
electromagnetic charge, and to the inverse mass, we may conclude that the contribution of very 
heavy quarks like, top and bottom, can be neglected compared with the light quarks (i.e. up, down, 
strange) and also the charm quark. The light quark masses are considered to be the same. Hence, 
the proton structure function may be written as 

F 2 = Ff" + Ff + Ff + Ff (3.6) 

Each component of the structure function associated with different quark flavors, is proportional to 
the particular quark wavefunction of the corresponding flavor, 

Ff(x,Q 2 ) = 2-^— [ d 2 b [ d 2 r [ dz\^ fl (r,z;m f ,e f )\ 2 N(r,b,x fl ) (3.7) 
47ra em J J J 

In the case of the charm quark, its contribution to the total structure function, has the following 
form 

Q 2 



'cm 



Ff(x,Q 2 ) = 2^^- I d 2 b I d 2 r \ dz\V cc (r, z;m c , e c )\ 2 N(r,b, x cc ) (3.* 



where the photon wavefunction, which decays into the charm-anticharm pair is given by 

|^(r,z;m c ,e c )| 2 = ^~{[* 2 + (1 + z) 2 ]a 2 Kl{ar) + m 2 c K 2 (ar) 

+ AQ 2 z 2 {\- z) 2 K 2 {ar)} (3.9) 

where a 2 is given by Eq. ( |3.5|) , and the numerical factor | is the squared value of the charge 
quantum number of the charm, in units of e. It is important to stress, that every quark flavor has 
an independent value of Bjorken-x, which depends on the quark mass. This point is presented in 
the next section. 
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3.2 Redefinition of Bjorken x 

The important ingredient of our approach, is the redefinition of Bjorken- a;. The motivation for this 
comes from the saturated region, or the region which lies below the saturation scale (Q 2 < Q 2 S ). In 
this region, the transverse momentum of the partons is proportional to the saturation scale, and 
this leads to the new definition of Bjorken-x 

4m 2 = (q + k) 2 

Am 2 = -Q 2 + 2P-q + k 2 = -Q 2 + 2P ■ q - Q 2 



Q 2 + Q 2 + 4m 2 . 



The kinematics is shown in Fig. |3.1| (fc) We can easily see, that this new definition of Bjorken-a; 
allows a smooth transition from large values of Q 2 , to the region of low values of Q 2 . At large Q 2 , 
we get back the ordinary expression for Bjorken-x, and at low values of photon virtuality, the main 
contribution comes from the saturation scale. 



3.3 Description of Fit 

From the fit to the experimental data, we can find values of the appropriate parameters for our 
model. For this purpose, we use all the latest data for the proton structure function F 2 from different 
collaborations in the region x < 0.01 and 0.045 < Q 2 < 150 GeV 2 . The small x cut leads to an 
upper limit on Q 2 [8-16]. It was observed, that the set of data from [8] should be rescaled with the 
factor of 1.037, in order to satisfy the best fit. For the numerical evaluation, we have implemented 
the formulae Eqs. (|3.2|) - fl3.4|) . We have also included the contribution from the heavy charm quark 
Eq. ( |3.8j ), and the new definition of Bjorken- x Eq. (|3.10| ). The model contains free parameters, 
whose values are to be determined from the fit to the experimental data. The result of our fit, with 
four quark flavors is represented in table 1. 





A 




^0 


Ql 


C 


m q 


m c 


X'ld.o.f. 


four quarks (u,d,s,c) 


0.785 


1.294 


0.060 


1.236 


1.0 


0.24 


1.3 


354/341 



Table 1: The parameters determined from the fit to the proton structure function F 2 



The parameters have the following physical meaning. A with ujq and C with /z 2 , and Qq, determine 
the hard scale and the initial gluon density, respectively. We observed, that the parameter C, has 
a strong correlation to other parameters. For every value of C, it is possible to find a set of four 
parameters, which satisfy almost the same x 2 - Hence, it was fixed with the arbitrary value C = 1.0. 
The value of the light quark mass, which corresponds to the best fit is m q = 0.24 (GeV). For the 
charm quark, the mass was taken to be m c = 1.3 (GeV). The reson for taking the same values for 
the three light quarks (u, d, s) originates from the fact, that the influence of the strange quark to 
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the total fit is relatively small, since its contribution in comparison with u and d is proportional to 
the ratio of charges squared. The resulting fit is plotted in Figs. [3.2| , |Q . 
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Figure 3.2: F2(x,Q 2 ) structure function of proton, as a function of Bjorken x for fixed value of photon 
virtuality Q 2 . This data was taken for small values of x < 0.01. These data points correspond to different 
collaborations. Asterisk corresponds to [10], diamond to [15] collaboration, pentagram to [13], triangles [8], 
circles correspond to [9] and hexagons to [14]- 
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Figure 3.3: (continued) F2(x,Q 2 ) structure function of proton, as a function of Bjorken x for fixed value 
of photon virtuality Q 2 . This data was taken for small values of x < 0.01 These data points correspond to 
different collaborations. Asterisk corresponds to [10], triangles [8], circles correspond to [9] and hexagons 

to [14]. 

4. Analysis of the model 

4.1 Comparison between the models 

The most widely spread saturation model, is based on the eikonal approach and has the following 
form 




(4.1) 



where Q is defined in Eq. ( |2.12| ). Using the AGK cutting rules [25], we can calculate the cross 
sections with different multiplicities k. This is the content of the AGK rules, and which relates 
the cross-section for observing a final state with /c-cut Pomerons, with the amplitudes for the 



exchange of m Pomerons F^> 



, oo 



T.^r-^T^-r^F^ (4.2) 



d?b ' kl(m — k)\ 

m=k 

We will apply this concept and calculate the predictions for the fc-cut cross-sections. We make a 
comparison between two models 



N 1 = and N n = 2 (l - e - *) (4.3) 



n 
+ 

For small values of Q, the dipole cross sections in Eq. ( [4.3|) are equal to Q, and are proportional 
to the gluon density. This allows one to identify the opacity, with the single Pomeron exchange 
amplitude of Fig. |4.1|. Hence, the multiple Pomeron amplitude is determined from the expansion 



i 




Figure 4.1: The single gluon-ladder contribution to the total 7* P cross section 

of the amplitudes in the form of a series expansion. For our model we obtain 

00 

N 1 = 2^(-l) m " 1 n m (4.4) 

m=l 

with 

F r (m) = Vl m (4.5) 
and for the eikonal model, the same procedure leads to 



where 



00 /o\ m 1 

*' = 2 5>ir- l g) if 

m=l x 7 



(4.6) 



*T=(* ~ { (4-7) 



The dipole cross section can be rewritten [26] in terms of F^ m \ as a sum over multi- Pomeron 



amplitudes 

CO 

N = 2^(-l) m - 1 F (m) (4.8) 

m=l 

The expression for the k cut Pomeron cross section, is obtained from the AGK cutting rules of 
Eqs.O, fl4j) and Q 



d 2 6 ^ v 7 fc!(m-A;)! l + 2Vt\l + 2Vt 



and 



M = V f _ ir ^ 2 ™ f^f — = — e" n (410) 

d 2 6 ^ l ; fc!(m-fc)! V 2 / m\ k\ V ' ' 

m=k ' 

The diffractive cross-section, is given by the difference between the total, and the sum over all 
cut cross sections 

d<Jdiff da tot _ sr^ dak (a-\i\ 

d?b ~ d?b ^ d?b [ ' ' 

fe=i 

and for two models reads as follows 

dfin = f SI \ _ JQ_ = 2tf 
d?b + l + Q (2fi + l)(ft + l) 1 ' ; 

and 

dal tf / \ / n\2 



= 2(l-e-§) -(l-e- n ) = (1-e-*) (4.13) 

Since we want to compare to different-valued functions, we need to normalize them. Finally, 
we compare between the ratios 

R * - l^Ft (4 ' 14) 

and 



where 



rM - da *l d% u 15) 

Hk ~ da*/d?b (4 ' 15) 



da P i v — -\ da u , . 

k=l 



13 



Below we present the plots describing the partial normalized cross sections (dak I d 2 b) / (da e i / d 2 b) , 
which correspond to fc-cut Pomerons and the diffractive cross section (dadiff/d 2 b). We compare the 
two models. Our approach, uses a solution of the non-linear evolution equation in a "toy-model" 
approach, and the Glauber approach. Our plots were calculated for fixed values of the dipole size 
r = 0.1/m, and the value Xf, = 0.001. Since higher cuts (k > 1) are strongly suppressed, they were 
rescaled by a factor 10 3fc . 



k = k=1 k = 2 k = 3 




b(fm) b(fm) b(fm) ' b (fm) 

Figure 4.2: Cross sections with different multiplicities as a function of impact parameter b (fm). Solid 
line corresponds to our model and dashed line to the eikonal approach. 



Although both models fit the DIS data quite well, there is a difference in the value, and in the 
shape of the partial cross sections, as shown in the plots. 

4.2 Shadowing corrections 

In this subsection, we want to discus the difference in the behavior of the two models (our model 
and the Glauber-like model), from the shadowing corrections point of view. Shadowing corrections 
(SC), are defined as the next to leading order terms in the series expansion of the amplitude near 
the amplitude zero point. 

N GF = = Slot ~ &o f + ~ - (4.17a) 



N Gl = l-e-% 1 = ^_i% + iSL_... (4.17b) 

2 4-2! 8-3! V ; 

We can easily conclude, that our proposed generating functional (GF) based model, and the 
Glauber-like (Gl) function, are the same at leading order, as was expected. Assuming this, we 
immediately obtain the next condition 

n Gt = 2tt Gf (4.18) 

Since we want to check the difference in the behavior of the two different parameterizations, we 
substitute the condition of Eq. (^4.18|) into Eq. (|4.17b|) . Finally, we obtain, the fact that our model, 



predicts larger SC than those proposed by the Glauber-like model. 

\ Gf = -Q 2 Gf + Q Gf 



SC GF = N GF -Q Gf =-Q 2 Gf + Q Gf -... (4.19a) 
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(4.19b) 



It is obvious that SC > SC . In spite of the fact that, generically, the SC are larger, we can 
fit all the experimental data, using the Glauber parametrization. It turns out that Qgu in such a 
parametrization, is larger than Qcf by 20 — 30%. The leading term is significant in the low energy 
domain, and should be consistent with the perturbative calculations of the amplitude. SC becomes 
significant at higher energies, especially in the saturation domain. Our model predicts the slower 
growth of the amplitude with energy, than that of the Glauber model. 



5. Predictions and descriptions within our model 

5.1 Charm quark contribution 

The developed dipole model, allows us to calculate a prediction for the inclusive charm quark 
contribution, to F 2 of the proton. Using the ansatz Eq. ( |3.8|) , we can easily obtain the contribution 
to the structure function from the charm quark F| c . 




Figure 5.1: Prediction of our model to F£ c ■ Solid line is our model, triangles correspond to HI 96-97 and 
circles to ZEUS 99-00 data [11, 12] respectively. Mass of the charm quark has the value: m c = 1.3 GeV 2 

Comparing this to the recent experimental data [11, 12] we find a good agreement between the 



calculated prediction, and the experimental result, as it can be viewed in Fig. [57L 

5.2 Description of dF 2 / d(\nQ 2 ) in HERA and the LHC kinematic region 

In this section, we want to check how our model is able to describe the slope of the structure 
function, as a function of the photon virtuality Q 2 , for fixed values of Bjorken-x, and vice versa. 
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For this purpose, we calculate the logarithmic derivative of F 2 

dF 2 

Since the only dependence on the virtuality Q 2 , is in the photon wavefunction, the resulting expres- 
sion for the calculation has the following form 



dF 2 



«9(lnQ 2 



Q 



2^2 

dQ 2 



d ( Q 



dQ 2 \A7r 2 a 



cr(x, Q 2 



F 2 + 



Q 4 da(x,Q 2 
4n 2 a dQ 2 



(5.2) 



and 



da(x, Q 2 
dQ 2 



2 I tfb I d 2 r I' dz d|*M^ 2 )l 2 Xt h) 
'o Jo Jo 



OO POO 

'2 l / j2„ 



(5.3) 



where \^(r, z, Q 2 )\ 2 and N(r, x, b) are defined in Eq. (|3.4|) and Eq. fl2.19|) , respectively. The resulting 
plots are presented in Figs. [5T2"| , |5T3"1 for fixed values of Q 2 and x respectively. Note, that these sets 
of experimental data, were not take into account in the fitting procedure. 
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Figure 5.2: Prediction of our model to logarithmic derivative dF2/d(lnQ 2 ) as a function of Bjorken-x for 
fixed values of Q 2 . 

We can see, that predictions which are based on our model, fits well with all experimental data 
on logarithmic derivative of F%. We enlarged the kinematic region towards very low x values, to 
give predictions for \q2 in the HERA and the LHC kinematic region. 
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Figure 5.3: Logarithmic derivative of F% as a function of Q 2 for fixed values of Bjorken-x. Solid line is 
our model and triangles correspond to HI 96-97 [8] respectively. 

5.3 Description of (91nF 2 /9(ln 1/x) 

In this subsection, we present our computation of = dlnF 2 /d(lnl/x). A comparison of our 
prediction with the HI experimental data [16] is shown in Figs. |5.4| , |5.5| for a fixed value of Q 2 and 
Bjorken-x respectively. 
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Figure 5.4: The logarithmic derivative X x = din i^/c^ln 1/x) as a function x for fixed values of photon 
virtuality Q 2 . Data taken from [16]. 
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Figure 5.5: The logarithmic derivative X x as a function of Q 2 for fixed values of Bjorken-x. 



We can see a good agrement with the experimental data. In order to investigate the behavior 
of the slope X x , in the high energy limit (very low x and Q 2 ), we expand our predictions to this 
region. The result is plotted in Fig. |5.6| . We want to pointed out the fact, that our model predicts 
a behavior of the structure function at the region of small photon virtualities Q 2 < 1 GeV 2 which 
is in good agreement with that obtained by Donnachie and Landshoff [27] from the fit to data. 
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Figure 5.6: The prediction for logarithmic derivative \ x = 9 In i^/c? (In 1/x) plotted at very low x and 

Q 2 . 



5.4 Prediction for F L at HERA 

Here, we want to present our predictions for the longitudinal part of the F 2 structure function F^. 
This longitudinal part, originates from the scattering of the longitudinally polarized virtual photon, 
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off the proton target. The expression for the calculation has the following form 



F L (x,Q 2 ) = ^2 J d 2 b j d 2 r J dz\^ L (r,z,Q 2 )\ 2 N(r,b,x) 



(5.4) 



where \^L(r,z,Q 2 )\ 2 corresponds to the longitudinal part of the photon wavefunction squared 
Eq. ( p.4|) , and N(r,b,x) is the interaction amplitude Eq. Q2.19| ). We used the relevant data on 
Fl from HI collaboration [8] in order to estimate our calculations. The experimental values of the 
longitudinal structure function, are not measured, rather they are extracted from the total structure 
function F 2 . The extraction of the longitudinal structure function, is based on the reduced cross 
section Eq. fl5.5|), which depends on two proton structure functions, F 2 (x, Q 2 ) and F L (x,Q2) 



a r = F 2 (x,Q 2 



[i + y) 



;F L (x,Q2) 



(5.5) 



where is y defined in Eq. ( |3.1|) . From the reconstruction of the kinematical variable, the desired 
data on Fl were obtained. Our main idea, is to predict the behavior of the longitudinal part of the 
structure function, in the HERA kinematic region (x < 10 -5 ), and to check how it fits the existing 
extracted data. The resulting plots are presented in Fig. |5T7| 
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Prediction for the Fl structure function at the HERA kinematics. The values of Q 2 are given 
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5.5 Diffractive production in DIS 

Diffractive deep inelastic scattering is usually characterized by two variables, the mass of the diffrac- 
tive system M%, and the momentum transfer t = (P — P') 2 . These variables are usually rewritten 
in terms of other, dimensionless variables as 



Q 2 + M 2 — t 



Q 2 + W 2 

which is the fractional energy-loss suffered by the incident proton and 

Q 2 



(5.6) 







Q 2 + M 2 -t 



(5.7) 



which corresponds to the momentum fraction, carried by a struck parton. The pomeron, which 
carries longitudinal momentum xp, is emitted by the proton, and subsequently undergoes hard 
scattering satisfying 



x B = x F (3 

The diffractive dissociation process is depicted in Fig. FT£ 
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Figure 5.8: Kinematic variables of semi-inclusive reaction ep — ► eNX diffractive dissociation 

The important property of the wavefunction formalism, is the ability to describe the diffractive 
scattering processes [28]. At small values of the diffractive mass M 2 , the elastic scattering of the qq 
pair, dominates, and the corresponding diffractive structure function reads as 



F 2 D (x,Q 2 ) T , L = I d% I rfV / l dz\* T ,L 

47r^a em J Jo 



[r,z;m f ,e f )\ N (r,b,x) 



(5.9) 



However, at larger values of the mass M 2 , the qqg contribution dominates, due to gluon pro- 
duction in the final diffractive state. We take into account the three leading twist terms 



,D(3) 



pi. + F L - + F 1 - 

qq 1 qq 1 qqg 



(5.10) 
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We introduce the last term in Eq. ( |5.1(J| ), to describe high mass diffraction, and as a simple ap- 
proximation of the first "fan" diagram, in which the emission of large numbers of gluons is taken 
into account. The longitudinal part Fk , has no leading logarithm in Q 2 , and can be neglected. 
The Feynman diagram for the interaction of a quark-antiquark dipole, with the target proton via 
two-gluon exchange, is shown in Fig. |5.9| . The gluons couple to the quark in all possible ways. 




Figure 5.9: qq contribution to diffractive DIS, where z represents the fraction of the energy of the photon 
that is carried by a quark. M is an invariant mass of diffractive system 



We follow the procedure proposed in Ref. [3,29] and we obtain 



QH4 /-oo r \M 2 -m 2 / 2 \ 3 



" ^IFK^ ~ HI JO JO . /I A 

and the transverse qq contribution has the form 

i ~T^ ) #(M,fla*0 + m 2 f <fi(k,b,l3,x F )} (5.12) 



where 



M 2 = Q2[l n ® with m 2 ± = k 2 ± + m 2 (5.13) 

P 



and the "impact factor" b, /3,x F ) (i = 0, 1, 2) given by: 



<j>i(k,b,p,xi>) = I drrKi \ \j k ^J f r | Ji(kr)N{x F ,r,b) (5.14) 

where K% and Jj are Bessel functions. This impact factor, represents the interaction between the 
produced dipole from the virtual photon, and the target. 

The next contribution qqg, was calculated assuming the strong ordering in the transverse mo- 
menta of the gluon and the qq dipole, namely k± g <C k± q> g. This assumption allows us to treat the 
qq, and qqg, as an effective color dipole in the transverse r space. The corresponding diagram is 
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plotted in Fig. [5.1U 



Pt 3 



Figure 5.10: gqg contribution with strong ordering with the gluon. Gluon's transverse momentum is 
considered to he smaller than the (anti)quark's one. 



Thus, we obtain: 

F^ g (i3,Q 2 ,xjp) 




W \- 2 a s (Q 2 ) [ x d0_ 
Z - 6 2vr J p (3' 



!QlT 2 X]p 
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0' 



+ 
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d 2 b I dk 2 In [ ^ 1 n P 1 ) <%(k,b,P',x P ) 











(5.15) 



Using the developed model (see section |2^2|) , which successfully fitted all the experimental data 
on DIS, we want to describe diffractive DIS, using the same model and parameters which were 
obtained from the fit. For this purpose, we used the latest data on diffractive dissociation, from 



ZEUS collaboration [14]. The resulting plots, are shown in Figs. |5.11|, |5.12 



6. Summary 

Using the generating functional approach, we develop a new saturation model, which describes well 
all the HERA data on deep inelastic scattering and diffractive dissociation. We took into account 
the QCD evolution, by evolving the initial gluon density with the DGLAP evolution equation. The 
innovation used in this current work, is the redefinition of Bjorken-x, in the saturation region, 
since the transverse momentum of partons cannot be neglected, and it is the saturation scale. We 
included in the calculation of the proton structure function the contribution from heavy quark 
production, and also investigated the behavior of the F| c structure function. The resulting x 2 /d.o.f 
for the fit of DIS data, is very close to one (x 2 / 'd.o.f = 354/341) and this reflects the fact that our 
new model is able to give reliable predictions. This was justified in the case of the description of 
diffractive dissociation data, using the same parameters which were obtained from the fit to DIS. 
Despite the fact that two completely different models are able to describe well the experimental 
data, we can distinguish between them by calculating the differential cross sections with various 
multiplicities do~ k /d 2 b. From figures shown, we can easily see the different behavior of these cross 
sections, as a function of the impact parameter b. From the plot of the saturation scale 27T , we can 
conclude, that up to the scale ~ 3 - 4 GeV 2 , saturation is significant, and the non-linear interaction 
term in the evolution equation, plays an essential role. We believe that this paper, introduces an 
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additional argument for the saturation phenomenon. It shows that the saturation models are able 
to describe all experimental data, including small values of Q 2 and low x. The description does not 
depend on the model assumptions widely used instead of the solution to the equation in the mean 
field approximation [30], which turns out to be rather complicated. 

7. Acknowledgments 

I would like to express my deep appreciation to Eugene Levin for his support in writing this paper. 
I am also very grateful to Asher Gotsman, Alex Prygarin, Jeremy Miller, and Alex Palatnik for 
fruitful discussions on the subject. Special thanks to Erez Etzion, Gideon Bella, Jonatan Ginzburg, 
Nir Amram and Eran Naftali for the technical support and essential discussions on the experimental 
background of the research. This research was supported in part by the Israel Science Foundation, 
founded by the Israeli Academy of Science and Humanities and by BSF grant # 20004019. 

References 

[1] A. H. Mueller, Nucl. Phys. B 415, 373 (1994). 

[2] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59 (1999) 014017 [arXiv:hep-ph/9807513]. 

[3] K. Golec-Biernat and M. Wusthoff, Eur. Phys. J. C 20 (2001) 313 [arXiv:hep-ph/0102093]. 

[4] J. Bartels, K. Golec-Biernat and H. Kowalski, Phys. Rev. D 66 (2002) 014001 
[arXiv:hep-ph/0203258]. 



[5] 


H. 


Kowalski and D. Teaney, Phys. Rev. D 68 (2003) 114005 [arXiv:hep-ph/0304189]. 


[6] 


K 


Golec-Biernat and S. Sapeta, Phys. Rev. D 74 (2006) 054032 [arXiv:hep-ph/0607276]. 


[7] 


H. 


Kowalski, L. Motyka and G. Watt, Phys. Rev. D 74 (2006) 074016 [arXiv:hep-ph/0606272]. 


[8] 


C. 


Adloff et al. [HI Collaboration], Eur. Phys. J. C 21, 33 (2001) [arXiv:hep-ex/0012053]. 


[9] 


S. 


Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 21 (2001) 443 [arXiv:hep-ex/0105090]. 


[10] 


J. 


Breitweg et al. [ZEUS Collaboration], Phys. Lett. B 487 (2000) 53 [arXiv:hep-ex/0005018]. 


[11] 


s. 


Chekanov et al. [ZEUS Collaboration], Phys. Rev. D 69 (2004) 012004 [arXiv:hep-ex/0308068] 


[12] 


c. 


Adloff et al. [HI Collaboration], Phys. Lett. B 528 (2002) 199 [arXiv:hep-ex/0108039]. 


[13] 


c. 


Adloff et al. [HI Collaboration], Nucl. Phys. B 497 (1997) 3 [arXiv:hep-ex/9703012]. 


[14] 


s. 


Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B 713 (2005) 3 [arXiv:hep-ex/0501060]. 


[15] 


M. R. Adams et al. [E665 Collaboration], Phys. Rev. D 54 (1996) 3006. 


[16] 


C. Adloff et al. [HI Collaboration], Phys. Lett. B 520 (2001) 183 [arXiv:hep-ex/0108035]. 



- 23 - 



[17] A. H. Mueller, Nucl. Phys. B 415 (1994) 373; ibid B 437 (1995) 107. E. Levin, Phys. Rev. D 49, 
4469 (1994). 

[18] E. Levin and M. Lublinsky, Nucl. Phys. A 730 (2004) 191 [arXiv:hep-ph/0308279]. E. Levin and 
M. Lublinsky, Phys. Rev. Lett. B 607 131 (2005) arXiv:hep-ph/0411121. 

[19] A. H. Mueller, Nucl. Phys. B 437 (1995) 107. 

[20] Yu. Kovchegov, Phys. Rev. D 60 (2000) 034008. 

[21] E. Levin and A. Prygarin, arXiv:hep-ph/0701178. 

[22] R. K. Ellis, Z. Kunszt and E. M. Levin, Nucl. Phys. B 420 (1994) 517 [Erratum-ibid. B 433 (1995) 
498]. 

[23] A. H. Mueller, Nucl. Phys. B 335 (1990) 115. 

[24] J.R. Forshaw and D.A. Ross, QCD and the Pomeron, Cambridge University Press, 1997; 
N. Nikolaev and B. G. Zakharov, Z. Phys. C 53 (1992) 331; 

V. Barone, M. Genovese, N. N. Nikolaev, E. Predazzi and B. G. Zakharov, Phys. Lett. B 326 (1994) 
161 [arXiv:hep-ph/9307248]. 

[25] V. A. Abramovsky, V. N. Gribov and O. V. Kancheli, Yad. Fiz. 18 (1973) 595 [Sov. J. Nucl. Phys. 
18 (1974) 308]. 

[26] A. H. Mueller and G. P. Salam, Nucl. Phys. B 475 (1996) 293 [arXiv:hep-ph/9605302]. 
[27] A. Donnachie and P. V. Landshoff, Nucl. Phys. B 244 (1984) 322. 
[28] N. Nikolaev and B. G. Zakharov, Z. Phys. C 53 (1992) 331. 

[29] E. M. Levin, A. D. Martin, M. G. Ryskin and T. Teubner, Z. Phys. C 74 (1997) 671 
[arXiv:hep-ph /9606443] . 

[30] M. Lublinsky, Eur. Phys. J. C 21 (2001) 513 [arXiv:hep-ph/0106112]; 

E. Levin and M. Lublinsky, Nucl. Phys. A 712 (2002) 95 [arXiv:hep-ph/0207374]; 
N. Armesto and M. A. Braun, Eur. Phys. J. C 20 (2001) 517 [arXiv:hep-ph/0104038]; 
K. J. Golec-Biernat, L. Motyka and A. M. Stasto, Phys. Rev. D 65 (2002) 074037 
[arXiv:hep-ph/0110325]; 

E. Gotsman, E. Levin, M. Lublinsky and U. Maor, Eur. Phys. J C 27 (2003) 411 
[arXiv:hep-ph/0209074]; 

E. Gotsman, M. Kozlov, E. Levin, U. Maor and E. Naftali, Nucl. Phys. A 742 (2004) 55 
[arXiv:hep-ph/0401021]; 

E. Iancu, A. H. Mueller and S. Munier, Phys. Lett. B 606 (2005) 342 [arXiv:hep-ph/0410018]; 



- 24 - 



1.2 GeV 



3 GeV 



to 




M x = 6 GeV 








11 GeV M x = 20GeV 



30 GeV 




0^- 
10 



10" 

X 




10" 
X 






10" 
X 



10 
X 




10" 
X 






Beta =0.0066 




10" 
X 



Figure 5.11: Diffractive dissociation structure function F 2 ((3,X]p,Q 2 ) as a function ofxp. Comparison of the prediction of the model 
with ZEUS data [14J. 




Figure 5.12: (continued) Diffractive dissociation structure function F 2 (/?,2?iP,Q 2 ) as a function ofxp. Comparison of the prediction 
of the model with ZEUS data [14]- 



